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Abstract 

We study the exit time r = T(o,oo) foi^ 1-dimensional strictly sta- 
ble processes and express its Laplace transform at as the Laplace 
transform of a positive random variable with explicit density. Conse- 
quently, r satisfies some multiplicative convolution relations. For some 
stable processes, e.g. for the symmetric |-stable process, explicit for- 
mulas for the Laplace transform and the density of r are obtained as 
an application. 



1 Introduction 

Let a G (0, 2) and (Xj, P^) be a strictly a-stable process in R with charac- 
teristic function 



[1 — i/3 arctan — sgn z 



E e""'^- = exp 

where /3 G [—1, 1] and /3 = for a = 1. For any C M let 



Td = inf{t >0: Xt^D} 

be the first exit time from D of the process Xf. Throughout this article we 
shall consider the starting point x > and 



T — T(0,oo), 
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the exit time of Xt from the positive half-hne. 

The question of the first exit time from domains are basic for all stochastic 
processes. Surprisingly few exact formulas are known for stable processes. 
The only exceptions are Brownian motion, completely asymmetric stable 
processes with a > 1 (see [H], [17]) and symmetric Cauchy process ([5], 
see also [H]). The quotient f/r was studied for independent r and dual f in 

Some recent results on this problem in the completely asymmetric case 
were obtained by T. Simon in |17i| and next were applied in [15] and [18]. On 
the other hand, some formulas were found by A. Kuznetsov in [12], however 
the final expressions are complicated. M. Kwasnicki in [T3] gives an integral 
representation of the density of r in the case of symmetric stable processes 
(/3 = 0). 

In this article we study the exit time r = r(o,oo) for 1-dimensional stable 
processes and give in Theorem [3] a new formula for its Laplace transform. 
It follows that r satisfies some multiplicative convolution relations (Corollary 
E]); in particular for a > 1 the exit time r is the multiplicative convolution of 
a l/a-stable subordinator with an explicitly given random variable p. We 
generalize in this way the result of [17j for all stable processes. Applications 
of Theorem [3] are next given in the final part of the article. New explicit 
formulas for the Laplace transform and the density of r are proven for the 
processes dual to those of the Doney's class Ci^i, in particular for the sym- 
metric |-stable process (Proposition [7] and Corollary [8]). Further applications 
of Theorem [3] will be presented in a forthcoming paper. 

The main tool to prove the results of this article is a series representation 
that we obtained in [9J for the logarithm of the bivariate Laplace exponent 
fi;(?7, 6) of the ascending ladder process built from the process Xf. This appli- 
cation of the series representation of In/t was announced in [H]. It allows to 
determine explicitely in Proposition [1] the inverse Stieltjes transform of the 
function 1/k(1, 6). 

2 Stieltjes transform and Wiener-Hopf fac- 
tors 

In this part of the article we will exploit our series representation of «;(!, 9) 
from [9] in inverting a Stieltjes transform. 
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Recall that if /i is a positive Borel measure on [0, oo) then for any x G 
(0, oo) the Stieltjes transform of /i is defined by 



whenever the integral converges. According to [3|, a function G on (0, oo) is 
of the form G{6) = a + Sfi{6) for a positive measure fi and a > if and only 
if 

(51) G extends to a holomorphic function in the cut plane C \ M_ 

(52) G{e) > for ^ > 

(53) lmG{z) < for Imz > 0. 
Then the inverse Stieltjes transform is 

S-\G){x) = -(1/vr) lim lmG{-x + iy) , x>0, 

where the limit, in general, is in vague sense and equals /i. If /i is absolutely 
continuous with a continuous density, the limit is equal to the density of fi 
for all X > 



Let a G (0,2) and {Xt,F^) be a strictly a-stable process in R. By 
i^a.piv^ ^) W6 denote the bivariate Laplace exponent of the ascending lad- 
der process built from Xt. We normalize it requiring that ^^,^(1,0) = 1. To 
simplify the notation we will write k{t],6) for a fixed pair a,p (or equiva- 
lently a fixed process Xt). By k we denote the Laplace exponent for the dual 
process Xt = —Xt. As usually we write the positivity coefficient 



p = P°(Xi > 0) = ^ + — arctan f/3tan ^] . 

2 na V 2 / 



Proposition 1. For p G (0, 1] \ {!/«} we have 
sin(pa7r) 1 x°'k{1,x) 



TT 



X + 6* + 2a;" cos(pa7r) + 1 k{1 



dx = (2) 



Proof. Denote G{6) = 1/k{1,6). The function G{6) extends to a holomorphic 
function hi{z) on C \ ]R_ (see p],(i) p. 205). Let C be the set of Liouville 
numbers. For 6 G (0, 1) and a ^ £ U Q we have by [9J 

/ ^ (-l)"+i^'"sin(pm7r) ^ (-1)'=+!^'= sin(paA;7r) 

G{e) = exp - 2. 2^ 



msin(— ) ^ ksiniakTc) 

m=l ^ a I fc=i V I 
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The right hand side of the last formula may be extended to a holomor- 
phic function h2{z) on {z G C: \z\ < 1} \ R_ defining = exp(aLogw) 
where Logw = ln|w| + iArgw, Aigw G {—11,71], is the principal value of 
the complex logarithm. We note that hi = h2 on (0, 1), hence hi = /i2 on 
{zGC: |z|<1} \ R_ and h2 extends to a holomorphic function on C \ ]R_, 
equal for \z\ > 1 to the holomorphic extension of ^^^-^y for > 1. 

In the first part of the proof we will compute 

l{x) = lim lmG{—x + iy) 

for positive x. Denote by h{z) the expression under exponential of h2. Let 
us compute for < x < 1 

1 1 
l{x) = lim lmexp{h{—x + iy)) = e^'^'-"'^ sin(Im(w)), 

TV y^0+ Vr 

where 

°° ' ^^™+i^-x)'"sin(pm7r) ^ (-l)'=+ie*"^^a;"'= sin(paA;7r) 



w 



msin(— ) ^ ksmiakn) 



The last limit is justified by a standard estimation argument, that implies 
that in a converging power series one can enter the limit under the series. 
Moreover, the same argument shows that when < a; < 1, we have 

l{x) = -- lim ImGH = --e^'=("')sin(ImH). (3) 

vr ui— ^— x,lm«)>0 TT 

Now we evaluate 

a;™sin(pm7r) (— 1)*^+^ cos(a/i;7r)a;"^ sin(paA;7r) 



^ sin pmTT v;^ -1 



msin(— ) ^-^ ksiniakn) 
J ^ ^ ^ (— 1)'^"'"^ sin(Q;A;7r)x"^ sin(pa;A;7r) ^ (— sin(pafc7r) 

fc=l ^ ' k=l 



We will need the following formulas from [TO 

m(kip) psincp 

= arctan 

k 1 — p cos ip 



y ^^^^ = arctan ^ G (0, 2vr), / < 1. (4) 



fc=l 
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f:^^^^ = 4Ml-2pcosv. + P^), v'G(0,27r),p^<l. (5) 

k=l 

Therefore applying a formula sin(arctan-u) = :^7j=^ we get 



• /T / w • / —x sm( pan) 

sinfim(w)) = sm(arctan -) 

1 + x° cos(pQ;7r) 



x" sin(pa7r) 

1+x" cos(pQ7r) _ -x" sin(pQ;7r) 



sin(pa7r) 



2 + 2x" cos(pa7r) + 1 ' 



cos(pQ7r) 

Now we compute Re(w). By (j5]) we get 

;-l)'"+ix™sm((l -p)m7r) 



Re(^) = ^ 



msin(— ' 

m=l ^ a ' 



^ sin((l - p)aA;7r) ^ {-l)''+^x'^^ cos{pak7i) 

< ^ ein I rvi"rr I < ^ 



/csinfafcvr) k 

k=\ ^ ' k=l 

log x) — - log(l + 2x" cos(pQ;7r) + x^"). 



Hence 



e^^^"") sin Im = — - — -— \ = ^ > 0. 

TT \ \ II ^ + 2x°cos(pa7r) + 1 ^ ' 

By [9l Lemma 5] we have for Q > \ 

and we use the same method and ([3]) to obtain 

/(x) = -l hm ImGH = ^^^^^ ^""^^f ^ , (6) 

TT to-s>-a;,lmw>0 TT X'^°' + 2x°' COS(pQ;7r) + 1 

for z > 0, X 7^ 1. 

As the function l{x) is continuous at x = 1 and by [8., p. 205] the limit 
lim^_>._i_im^t,>o ImG'(w) exists, it follows that the convergence in ([6]) holds also 
for X = 1. 
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Let us now justify the fact that the function G{9) = 1/ 9) is a Stieltjes 
transform of a positive measure on R"*". We will check the conditions (Sl-3) 
given in the beginning of this section. 

The function G{d) is strictly positive for 9 G (0, oo) and it extends to a 
holomorphic function hi{z) on C\M_. Thus the conditions (SI) and (S2) 
are verified. In order to justify (S3), we use the following property that we 
proved above: the harmonic function — (l/7r)ImG extends continuously to 
the closed upper half-space {Imz > 0} and its boundary values on R are 
l{—x) > when x < and for x > 0. Taking into account the fact that 
\im\z\-^oo G{z) = ([8, p. 205]), the maximum principle([2, 1.10]) implies that 
lmG{z) < on {Im^ > 0} and (S3) also holds. 

It follows that for a certain a > we have G{9) = a+S{l){9). Considering 
— 7- oo we determine a = 0. 

Finally consider any a G (0,2]. Since the Lebesgue measure of the set 
£ U Q is we can take a sequence a„ tending to a. Passing to the limit we 
obtain ([2]) for all a G (0,2]. □ 

Remark. Other proofs of the fact that 1/k(1, 9) is the Stieltjes transform 
of a positive measure /i seem possible, using properties of Bernstein functions 

m)- 

We deduce immediately from Proposition [T] the following corollary. 
Corollary 2. For p E [0, 1) \ {1 — 1/a} we have 
sin((l - p)a7r) p 1 x"/€(l,a;) 



dx = ^ . (7) 



TT 7o X + 61x2" + 2x"cos((l - p)a7r) + 1 k(1, 

3 Laplace transform of r and applications 

The following theorem is the main result of the article. 

Theorem 3. Let Xt be a non- spectrally positive strictly a-stable process on 
M. For any t > we have 

,1 sin((l - p)a7r) .i/^^ x'^-^k{1,x) 



vr Jo x^" + 2x" cos((l - p)a7r) + 1 

Remark. Observe that the only case excluded from the Theorem [3] 
is well known: when Xt is a spectrally positive a-stable process starting 
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from Xq = X, 1 < a < 2, then (y'T(Q^))x>o is a 1/a-stable subordinator and 
]gig-iT _ g-i^''" Q4j p.281). When Xt is spectrally negative, the formula ([8]) 
was obtained recently by T.Simon([17J). 

Proof. We note that if ^'(0) = 1 then 



oo 



-/(x)c/x = g{t) =^ / :!^dx = -- (9) 



Indeed 



Jq x + t X Jq x + t X 

From [H] we know that 

where r = T(o,oo) and r],9 > 0. Putting rj = 1 and applying ([9]) to ([7]) we get 
sin((l - p)a7r) p 1 



(10) 



TT 7o a; + 6'x2° + 2x'^cos((l-p)a7r) + l 



sin((l - p)a7r) p r .-v{x+e) a:" ^^^^ 



oo 



TT 7o Jo + cos((l - p)a7r) + 1 



vr 7o + 2x" cos((l - p)a7r) + 1 ) 



Therefore 



TT 7o + 2a;" cos((l - p)a7r) + 1 

and assertion of the theorem follows from scaling property E^e^^ = E^e~^"^ 
of stable processes. □ 

For a > 1 we immediately obtain from Theorem [3] a formula for the 
density h = ha,p of r under P-*^. 
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Corollary 4. Let a G (1,2). The density of r = T(o,oo) under P" is given by 

sm((l - p)a7r) f°° x''-^k{1,x) 

his) = / rii/a[x,s)— ^ ^ ax, 11 

^' vr Jo + 2x°cos((l -p)a7r) + 1 ' ^ ' 

where rj^{t,x) is the transition density of j-stable suhordinator. 

Proof. Since e~^^ri^{t,x)dx = e~^^'' we obtain by Theorem [3] 



St 



e ^^hJs) ds 







sin((l-p)a7r) x-i«:(l,a:) , , 
/ / e ^*r/i/,(x,g) „ , ———dsdx 



71" Jo JQ ^ 

oo 



+ cos((l - p)a-n) + 1 



sin((l - p)a7r) [°° x''-^k{1,x) 



' e / ?7i/„(x, s)— r — — ; : dx ds 

TT Jq ' x2" + 2a;"cos((l-p)a7r) + l 

and the assertion follows. □ 

Corollary 5. If Xt is a Cauchy process on M (a = l,p = 1/2^) then the 
density of the exit time r = T(o,oo) under is given by 

, , , 1 x) 

IT X^ + 1 

Remark. The above formula for /ii,i/2 was obtained previously by Darling 
in [5j (see also [m (7.13)]). 

For a 7^ 1, Theorem [3] gives interesting multiplicative convolution relations 
verified by r. We present them in the following subsection. 

3.1 Interpretation in terms of multiplicative convolu- 
tions 

For a given strictly a-stable process Xf with p G [0, 1) \ {1 — l/a} we define 
the following function on M"*" 

_ sin((l - p)a7c) ft:(l,x^/") 

"""•"^^^ Tra x2 + 2x cos((l - p)a7r) + 1 ' ^ " 
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Observe that ma^p{x) is a probabihty density on M"*". This follows from the 
formula 

sin((l - p)«7r) 

ax = 1 



TT Jq + 2x"cos((1 - p)a7r) + 1 

obtained from ([8]) when t — )■ and by a change of variables x = y^^°'. 

Denote by Ma,p a positive random variable with density rria^p. The vari- 
able Mq, i/a, 1 < a < 2 appeared for the first time in [17J for the special case 
of a completely asymmetric a-stable process, 1 < a < 2 (in our context a 
spectrally negative process), when k(1,x^/") = 1 + x^^°'. 

Let 7 G (0,1) and ri^(t,x) be the transition density of 7-stable subor- 
dinator. Denote by A^(7) a random variable with the density ri^{l,x), i.e. 
Eexp(-a;A^(7)) = 6"^^ x > 0. 

Recall that if Y and Z are independent random variables on (0, 00) with 
densities / and g respectively, then the multiplicative convolution Y x is 
a random variable with the density 

nYxZ^edt]= r/(^)^?M§- (12) 

Corollary 6. (i) Let 1 < a < 2. Suppose that random variables Ma^p and 
N{l/a) are independent and that Xq = 1. We have 

T ^ M^^p X N{l/a), 

(a) Let < a < 1. Suppose that random variables t and N{a) are indepen- 
dent. We have 

T X N{a)'' = M^^p. 

Proof. Part (i) follows immediately from Corollary HI In f fTTj) we substitute 
= M and use the scaling property rji/aiu^^", s) = u~^T]i/a{^, su~^). 
In order to prove (ii), we use e~*^ = E,[—(tyy^°'N{a)] in the left-hand side 

of ([H]), we apply Fubini and change variables x = y^^^'u. By unicity of the 

Laplace transform we get 

x^ du sin((l — p)a7r) K(l,a;) 

h^[—)ri^[u) 



'u"" u"' na a;2" + 2x"cos((l - p)a7r) + 1 

Replacing by x in the last formula and using f fT2|) ends the proof of (ii). □ 
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3.2 Application for the Doney's class Ci i 

Let a G [|, 1) and 1 — p = l/a — 1, i.e. we consider a process dual to a 
process from the class C\^\ from R. Doney's article [7]. We denote the class 
of such processes by C\^\. In this case we have by [7] or [9] 

x^" — 2a;" cos avr + 1 
.(l,x) = . 

We recall the definition of the confluent hypergeometric functions \F\ and 
U (see [1]) 

fc=0 ^ ^'^ 

[/ (a, 6, 2;) 

TT 



iFi(a,6,2) + g - &,2 - 6,2:) ;2 G C \ 



sinvrfe Vr(l + a-6)r(6) r(a)r(2 - 6) 

where (a)^ = a(a + 1) ... (a + A; — 1), (a)o = 1 is Pochhammer symbol. Using 
Theorem [3] we get the following formulas for Laplace and Stieltjes transforms 
of r 

Proposition 7. If a E [i, 1) and Xt G Ci^i then 

(,) Ele--* = ^^^^^(a)^(l-a,^l/°)e*''^ (13) 
vr 

^1 1 r^-. A--^"-"^" .V"-SF,(l,2-c.,.V"x-V") \ 

(14) 

where V[a,z) = t"'~^e~* dt is the incomplete Gamma function. 
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Proof. To prove (i) we use simple transformations of integrals 
,1 sin(a(l/a - l)7r) /"^ e-*'''"^^"-^^^" - 2a;" cos avr + 1) 



E^e" = ^ ^ ' — / ^ ■ - dx 

n Jo (l + x)(x2" + 2x°cos(a(l/a-l)7r) + l) 



TT Jo 1 + X 



(ix (15) 



sm avr 



TT 




^0 



sin(a7r) p r(a) ^ sin(a7r)r(a) /"^ ^^i/..,^ 



= ^^r(a)r(l-a,tV")e*^^". 

TT 

To prove (iz) we use the following integral representation of U (see [1]) 

1 /"°° 

U(a,b,z) = —- e-'Y-^l + tf-^-Ut, Rez>0. 
r(a) Jo 

Applying this to ( IT5|) we get 



x + r Jo Jo 1 + s 



Jo Jo a;(l + s) 

sinfavr 



e r(a) du. 



TT Jo X 

Hence 



^^ = 70 — fi2^) 



_„ , e- - ni/"-SFi(l,2 - a,Mi/"a;-i/") \ , 
e 71 TT du. 

□ 

It is possible to invert Stieltjes transform in ( fT4l) and the resulting density 
of T is given by a series either in x or in 1/x (cf. [12j). For a = 2/3 the 
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process Xt is symmetric and the density of r has a nice integral representation 
involving hypergeometric function 

Corollary 8. Let Xt he a symmetric "^-stable process on M with Xq = 1. 
Then the density of r = T(o,oo) is given by the formula 

and its Laplace transform is 

E^e-*^ = ^r(2/3)r(l/3, t'/^)e'''' , t > 0. 
Proof. The second part is a direct consequence of (fT3l) . By ( fT4l) 

^ x + T~ J, [ X x3/2 r(4/3) ' 
Inverting Stieltjes transform we get 

F\t G dx) 

- _i im ^ e"''"^"''"^"''^^' u'/^ 4/3, M3/2^-3/2e-3W2) 



TT 



-X a;3/2e3W2 r(4/3) 

1 l^ e-^^^--^^^- ^V SF,(1, 4/3, ^3/2^-3/^,, 

TT io I a;'/' r(4/3) ' 



Since 



Im2iFi(l,4/3,?/i) = Re^ 



k=0 



fc!(4/3)fc 



k=0 



,2\k 



(4/3) 



2fc 



E 

fc=0 



,2\A; 



fc!(4/6)fc(7/6),22^ 



LF2(l;2/3,7/6;-y74) 
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we finally get 

P^(r G dx) 



1 p / sin(n3/2x'3/2) ni/2^F2(l;2/3,7/6;-n3a:-V4) ' 
V ^ ^ r(4/3) ' 



TT 







1 r,..(,„(,3/.)^ ^-^^.(l;^2A7/6;-.V4) ),^. 



TT 



□ 

Remark. It is possible to obtain the formula ( !T6|) from the results of 
M. Kwasnicki [13]. However, in order to do this, one has to make several 
non-elementary transformations of integrals and our approach seems simpler 
than using 
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